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Superconductors in realistic geometries: Geometric edge barrier versus pinning 
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The magnetic response of type-II superconductors can be irreversible due to two different rea- 
sons: vortex pinning and barriers for flux penetration. Even without bulk pinning and in absence 
of a microscopic Bean-Lingston surface barrier for vortex penetration, superconductors of nonellip- 
soidal shape can exhibit a large geometric barrier for flux penetration. This edge barrier and the 
resulting irreversible magnetization loops and flux-density profiles are computed from continuum 
electrodynamics for superconductor strips and disks with constant thickness, both without and with 
bulk pinning. Expressions for the field of first flux entry H en and for the reversibility field H ICV 
above which the pin-free magnetization becomes reversible are given. Both fields are proportional 
to the lower critical field H c i and else depend only on the specimen shape. These realistic results 
are compared with the reversible magnetic behavior of ideal superconductor ellipsoids. 

PACS numbers: 74.60. Ec, 74.60.Ge, 74.55.+h 
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I. INTRODUCTION 



The irreversible magnetic behavior of type-II super- 
conductors usually is caused by pinning of the vortices 
at inhomogeneities in the material (!]]. However, similar 
hysteresis effects were also observed in type-I super- 
conductors, which do not contain flux lines, and in type- 
II superconductors with negligible pinning. In these two 
cases the magnetic irreversibility is caused by a geomet- 
ric (specimen-shape dependent) barrier which delays the 
penetration of magnetic flux but not its exit. In this 
respect the macroscopic geometric barrier behaves sim- 
ilar as the microscopic Bean-Livingston barrier |j| for 
straight vortices penetrating at a parallel surface. In both 
cases the magnetic irreversibility is caused by the asym- 
metry between flux penetration and exit. The geometric 
irreversibility is most pronounced for thin films of con- 
stant thickness in a perpendicular field. It is absent only 
when the superconductor is of exactly ellipsoidal shape or 
is tapered like a wedge with a sharp edge where flux can 
penetrate easily due to the large local enhancement of 
the external magnetic field at this edge in a diamagnetic 
material. 

Ellipsoids are a particular case. In superconducting 
ellipsoids the inward directed driving force exerted on 
the vortex ends by the surface screening currents is ex- 
actly compensated by the vortex line tension B. An 
isolated vortex line is thus in an indifferent equilibrium 
at any distance from the specimen center. The repulsive 
vortex interaction therefore yields a uniform flux density 
and the magnetization is reversible. However, in speci- 
mens with constant thickness (i.e. with rectangular cross- 
section) this line tension opposes the penetration of flux 
lines at the four corner lines, thus causing an edge bar- 
rier; but as soon as two penetrating vortex segments join 
at the equator they contract and are driven to the speci- 
men center by the surface currents, see Figs. 1 and 2. As 



opposed to this, when the specimen profile is tapered and 
has a sharp edge, the driving force of the screening cur- 
rents even in very weak applied field exceeds the restoring 
force of the line tension such that there is no edge bar- 
rier. The resulting absence of hysteresis in wedge-shaped 
samples was clearly shown by Morozov et al. [|J. 

For thin superconductor strips with an edge barrier an 
elegant analytical theory of the field and current profiles 
has been presented by Zeldov et al. using the theory 
of complex functions, see also the calculations • With 
increasing applied field H a , the magnetic flux does not 
penetrate until an entry field H en is reached; at H a = H en 
the flux immediately jumps to the center, from where it 
gradually fills the entire strip or disk. This behavior in 
increasing H a is similar to that of thin films with ar- 
tificially enhanced pinning near the edges [0J5), but in 
decreasing H a the behavior is different: In films with en- 
hanced edge pinning (critical current density J® dge ) the 
current density J at the edge immediately jumps from 
_l_jodgc £ Q _jodgc w \ ien th e ramp rate reverses its sign, 
while in pin-free films with geometric barrier the current 
density at the edge first stays constant or even increases 
and then gradually decreases and reaches zero at H a = 0. 
For pin-free thin strips the entry field H en was estimated 
in Refs. §,0Jll). 

The outline of the present work is as follows. Section 
2 discusses the reversible magnetic behavior of pin-free 
superconductor ellipsoids. The effective demagnetization 
factor of long strips (or slabs) and circular disks (or cylin- 
ders) with rectangular cross section 2a x 2b is given in 
Sec. 3. In Sec. 4 appropriate continuum equations and 
algorithms are presented that allow to compute the mag- 
netic irreversibility caused by pinning and/or by the ge- 
ometric barrier in type-II superconductors of arbitrary 
shape, in particular of strips and disks with finite thick- 
ness. Results for thick long strips and disks or cylinders 
with arbitrary aspect ratio b/a are given in Sec. 5 for 
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FIG. 1. Field lines of the induction B(a;, y) in strips with 
aspect ratio b/a = 2 (top) and b/a — 0.3 (bottom) in perpen- 
dicular magnetic field H a . Top left: H a /H c \ = 0.66, in in- 
creasing field shortly before the entry field H cn /H c i — 0.665. 
Top right: H a /H c i = 0.5, decreasing field. Bottom: 
Ha/H c \ = 0.34 in increasing field just above H cn /H c i = 0.32. 
Note the nearly straight field lines in the corners indicating 
the tension of the flux lines. The field lines of cylinders look 
very similar. 

pin-free superconductors and in Sec. 6 for superconduc- 
tors with arbitrary bulk pinning. In particular, explicit 
expressions are given for the field of first flux entry H en 
and for the reversibility field H TCV above which the mag- 
netization curve is reversible and coincides with that of 
an ellipsoid. 



II. ELLIPSOIDS 

First consider the known magnetization of ideal 
ellipsoids. If the superconductor is homogeneous 
and isotropic, the magnetization curves of ellipsoids 
M(H a ; N) are reversible and may be characterized by 
a demagnetizing factor N. If H a is along one of the 
three principal axes of the ellipsoid then N is a scalar 
with < N < 1. One has N — for long speci- 
mens in parallel field, N = 1 for thin films in perpen- 
dicular field, N = 1/2 for transverse circular cylinders, 
and N — 1/3 for spheres. For general ellipsoids with 
semi-axes a, b, c along the cartesian axes x, y, z, the 
three demagnetizing factors along the principal axes sat- 
isfy N x + N y + N z = 1. 
a = b one has 7V r = N,, 



with a — b < c and for disks with a = b > c with eccen- 
tricity e=\l — c 2 /a 2 ) 1 / 2 one obtains (l2) 



1 



N a = 



-(atanhe — e), (cigar), 
-(e — atane), (disk). 



For thin ellipsoidal disks with a > b^$> c one has |1J 

c . 



1 



-E(k), 



(1) 



(2) 



where E(k) is the complete elliptic integral of the second 
kind with k 2 = 1 — b 2 /a 2 . 

When the magnetization curve in parallel field is 
known, M(H a ; 0) = B / fio — H a where B is the flux den- 
sity inside the ellipsoid, then the homogeneous magneti- 
zation of the general ellipsoid, M(H a ; N), follows from 
the implicit equation 



Hi = H a - N M(Hi] 0) . 



(3) 




For rotational ellipsoids with 
= (1 — N z )/2 where for cigars 



FIG. 2. 3D plots of the screening current density J 3 (x,y), 
Eq. (11), in superconductor strips with b/a — 2 (top) and 
b/a = 0.3 (bottom) as in Fig. 1. Shown is the limit of small 
applied field H a <^ H c \ before magnetic flux has penetrated. 
For better presentation the depicted J 3 (x,y) is smeared over 
a few grid cells. 
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Solving Eq. (3) for the effective internal field Hi, one ob- 
tains M = M(H a ;N) = M(Hi;0). In particular, for the 
Meissner state (B = 0) one finds M(H a ; 0) = — H a and 



M(H a ; N) 



1 — N 



for \H a \ < (l-N)H cl . (4) 



At the lower critical field H c i one has Hi = H cl , H a — 
H' cl = (1 - N)H cl , B = 0, and M = -H cX . Near the 
upper critical field H C 2 one has an approximately linear 
M(H a ; 0) = j(H a - H c2 ) < with 7 > 0, yielding 



M(H a ;N) 



7 



1 +jN 



(H a -H c2 ) for H a ~H c[ 



(5) 



Thus, if the slope 7 <C 1 is small (and in general, if 
\M/H a \ <C 1 is small), demagnetization effects may be 
disregarded and one has M(H a ; N) w M(H a ; 0). 

The ideal magnetization curve of type-II supercon- 
ductors with N = 0, M(H a ;0) or B(H a ;0)/^ o = 
H a +M (H a ; 0), may be calculated from Ginzburg-Landau 
(GL) theory |14|, but to illustrate the geometric barrier 
any other model curve may be used provided M(H a ; 0) = 
—M(—H a ; 0) has a vertical slope at H a = H c \ and de- 
creases monotonically in size for H a > H c \. Below for 
simplicity I shall assume H c \ <C H c2 (i.e. large GL pa- 
rameter k 3> 1) and H a <C H c2 . In this case one may use 
the model M(H a ;0) = -H a for \H a \ < H cl and 



M{H a ;0) = {H a /\H a \)(\H a \ A - H 6 cl ) 



3 Nl/3 



H a 



(G) 



for \H a \ > H c i, which well approximates the pin-free GL 
magnetization Jljj ]. 



III. THICK STRIPS AND DISKS IN THE 
MEISSNER STATE 

In noncllipsoidal superconductors the induction B(r) 
in general is not homogeneous, and so the concept of a 
demagnetizing factor does not work. However, when the 
magnetic moment m = i J r x J(r)d 3 r is directed along 
H a , one may define an effective demagnetizing factor N 
which in the Meissner state (B = 0) yields the same slope 
Mj H a = -1/(1-N), Eq. (2), as an ellipsoid with this AT. 
Here the definition M — m/V with m = mH a /_ff a and 
specimen volume V is used. In particular, for long strips 
or slabs and circular disks or cylinders with rectangular 
cross-section 2a x 2b in a perpendicular or axial magnetic 
field along the thickness 2b, approximate expressions for 
the slopes M/H a = m/(VH a ) are given in Refs. @Jl|]. 
Using this and defining q = (\M/H a \ — l)(b/a), one ob- 
tains the effective N for any aspect ratio b/a in the form 



N = 1- 1/(1 + qa/b) 



^strip , 

4 

3^ 



0.64 tanh 
2 



<Zdisk 



3tt 



tanh 



0.64- In (1.7 
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1.27- In ( 1 
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FIG. 3. Top: The axial magnetic induction B y (r,y) in the 
midplane y — of a pin-free superconductor disk with aspect 
ratio b/a = 0.3 in increasing field (solid lines) and then de- 
creasing field (dashed lines), plotted at H a /H a \ = 0.4, 0.42, 
0.5, 0.52, 0.6, 0.7, 0.8, 0.7, 0.6, 0.1, 0. Bottom: The 
induction B y (0,0) in the center of the same disk (solid line) 
and of a strip (dashed line), both with b/a = 0.3. The sym- 
bols mark the field values at which the profiles are taken. 
Also shown is the magnetization loop for the same disk and 
strip and the corresponding reversible magnetization (dotted 
lines) . 

In the limits b <C a and b a, these formulae are exact, 
and for general b/a the relative error is < 1%. For a = b 
(square cross-section) they yield for the strip N = 0.538 
(while N = 1/2 for a circular cylinder in perpendicu- 
lar field) and for the short cylinder N — 0.365 (while 
N = 1/3 for the sphere) . 



IV. COMPUTATIONAL METHOD 



To obtain the full, irreversible magnetization curves 
M(H a ) of nonellipsoidal superconductors one has to re- 
sort to numerics. Appropriate continuum equations and 
algorithms have been proposed recently by Labusch and 
Doyle 1 17 1 and by the author Eq], based on the Maxwell 
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FIG. 4. The irreversible magnetization curves —M(H a ) of 
pin-free circular disks and cylinders with aspect ratios b/a = 
0.08, 0.15, 0.25, 0.5, 1, 2, 5, and oo in an axial field (solid 
lines). Here the irreversibility is due only to a geometric 
edge barrier for flux penetration. The reversible magneti- 
zation curves of the corresponding ellipsoids defined by Eqs. 
(3,6,7) are shown as dashed lines. 

equations and on constitutive laws which describe flux 
flow and pinning or thermal depinning, and the equilib- 
rium magnetization in absence of pinning, M(H a ; 0). For 
arbitrary specimen shape these two methods proceed as 
follows. 

While method [jlTj considers a magnetic charge den- 
sity on the specimen surface which causes an effective 
field Hi(r) inside the superconductor, our method |l8| 
couples the arbitrarily shaped superconductor to the ex- 
ternal field B(r,i) via surface screening currents: In a 
first step the vector potential A(r, t) is calculated for 
given current density J; then this linear relation (a ma- 
trix) is inverted to obtain J for given A and given H Q ; 
next the induction law is used to obtain the electric field 
[in our symmetric geometry one has E(J, B) = —dA/dt], 
and finally the constitutive law E = E(J,B) is used to 
eliminate A and E and obtain one single integral equa- 
tion for J(r,i) as a function of H a (t), without having 
to compute B(r, t) outside the specimen. This method 
in general is fast and elegant; but so far the algorithm 
is restricted to aspect ratios 0.03 < b/a < 30, and to a 
number of grid points not exceeding 1400 (on a Personal 
Computer). Improved accuracy is expected by combin- 
ing the methods [jlTj (working best for small b/a) and 
]i~8[ . Here I shall use the method jl8| and simplify it to 
the two-dimensional (2D) geometry of thick strips and 
disks. 

In the 2D geometry of thick strips Jl5|] or short cylin- 
ders |l(| in an applied magnetic field B a = ^i H a = V x 
A a along y, one writes r = (x, y) or r = (p, y) (in cylin- 
drical coordinates p, (p, y) . For a homogeneous applied 
field the applied vector potential in these two geome- 
tries reads A a — —xB a or A a = —pB a /2. The current 



density J(r,i), electric field E(r, i), and vector potential 
A(r, t) now have only one component oriented along z or 
ip and denoted by J, E, A. The method HJjJIJ de- 
scribes the superconductor by its current density J(r,t), 
from which the magnetic field B(x,y,i) = (B x ,B y ) or 
B(p, y,t) = (Bp, By), the magnetic moment m(t) (along 
y), and the electric field E(r,t) = E( J, B,r') follow di- 
rectly or via the constitutive law E = E(J, B). For high 
inductions B ^> poH c i one has B = /irjH everywhere 
and J = —V 2 (A — A a ). The current density J is then 
obtained by time-integrating the following equation of 
motion, 



j(r,t) 



d 2 r' K(t,t')[E(J,B) + A a (v' ,t)} 



(8) 



Here K(r, r') = <5(r,r') _1 is an inverse integral kernel 
obtained by inverting a matrix, see |l5|,[l6| for details. 
The kernels Q and K apply to the appropriate geome- 
try and relate J to the current-caused vector potential 
A — A a in the (here trivial) gauge V • A = via integrals 
over the specimen volume V , 



A(r) = -f dVQ(r,r')J(r')U(r), 
Jv 



(9) 



J(r) 



dVf(r,r')[i(r')-i(r')]. (10) 



The Laplacian kernel Q is universal, e.g., Q(r,r') — 
— (/in/27r) In |r — r'| for long strips with arbitrary cross 
section, but the inverse kernel K depends on the shape 
of the specimen cross section. Putting A(r') — in Eq. 
(10) (Meissner state) one sees that 



J s (v)= [ d 2 r'K(r,r')A a (r') 
Jv 



(11) 



is the surface screening current caused by the applied 
field. In particular, one has J s (r) = inside the super- 
conductor. In our method J s automatically is restricted 
to the layer of grid points nearest to the surface, see Fig. 
2. 

If one is interested also in low inductions one has to 
generalize Eq. (8) to general reversible magnetization 
H = H(B). This is achieved by replacing in the con- 
stitutive law E(J, B) the genuine current density J = 

1 V x B by the effective current density Jh = V x H 
which drives the vortices and thereby generates an elec- 
tric field E. That Jh = V x H(B, r) enters the Lorentz 
force is rigorously proven by Labusch 17 1. Within the 
London theory this important relation may also be con- 
cluded from the facts that the force on a vortex is deter- 
mined by the local current density at the vortex center, 
while the energy density F of the vortex lattice is deter- 
mined by the magnetic field at the vortex centers. Thus, 
Jh = V x (dF/ <9B) is the average of the current densities 
at the vortex centers, which in general is different from 
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FIG. 5. The magnetization curves M(—H a ) = -M(H a ) of 
a thick disk with aspect ratio b/a = 0.25 and constant H c i 
for various pinning strengths, J c — 0, 0.25, 0.5, 1, 1.5, 2, 3, 4 
in units H c i/a, and various sweep amplitudes. Bean model. 
The inner loop belongs to the pin- free disk ( J c — 0), the outer 
loop to strongest pinning. The reversible magnetization curve 
of the corresponding ellipsoid is shown as a dashed curve. 

the current density J = /1q [ VxB averaged over the vor- 
tex cells. In our 2D geometry one thus has to replace in 
Eq. (8) 



E[J(r'),B(r')]^E[J H {r'),B(r')], 



(12) 



where Jh = dH y /dx — dH x /dy depends on the reversible 
material law H(B) = dF/dB with H x = H(B)B X /B, 
H y = H(B)B y /B, and B = (B 2 C + B*) 1 / 2 . 

The boundary condition on H(r) is simply that one 
has H = B//io at the surface (and in the vacuum out- 
side the superconductor, which does not enter our calcu- 
lation). This boundary condition may be forced by an 
appropriate space-dependent material law H = H(B,r), 
which outside and at the surface of the superconductor 
is trivially H = B//io- The specimen shape thus enters 
in two places: via the integral kernel K (r, r') and via the 
material law H = H(B, r). 

To compute the induction B(r) entering H(B), for 
maximum accuracy one should not use the derivative 
B = V x A but the Biot-Savart integral 



B(r)= / dVL(r,r')J(r')+B a (r) 



(13) 



with appropriate kernel L(r, r'). The accuracy of the 
method then depends mainly on the algorithm used to 
compute the derivative Jh = V x H. A useful trick 
is to compute Jh as Jh = J + V x (H — B/^o) where 
H — B//io is typically small and vanishes at the surface. 




FIG. 6. Magnetization curves of a disk as in Fig. 5 but with 
J c =const. and for various lower critical fields H c \ = 0, 0.2, 
0.35, 0.5, 0.7, 1 in units aj c . Bean model. 

For the following computations I use simple models for 
the constitutive laws of an isotropic homogeneous type-II 
superconductor without Hall effect, though our method 
|lj) is more general. With Eq. (6) and H = B/fi — M 
one has 



H(B)=^[Bl 1+ B^/\ 



(14) 



with B c i — /io-ffci- A simple independent current - 
voltage law which describes pinning, thermal depinning 
and flux flow is E(J, B) = p(J, B)3 with 

{J/JcY 



p(J, B) = p B 



i + (J/JcY 



(15) 



This model has the correct limits p oc J a (J <C J Cl flux 
creep) and p = pqB = ppp (J 3> J c , flux flow, po = 
const), and for large creep exponent a ^> 1 it reduces to 
the Bean critical state model. In general the critical cur- 
rent density J c = J C (B) and the creep exponent a(B) > 
will depend on B. For pin-free superconductors ( J c — > 0) 
this expression describes usual flux flow, i.e., viscous mo- 
tion of vortices, E = ppp(-B)J, with flux-flow resistivity 
/Off oc B as it should be. 



V. PIN-FREE SUPERCONDUCTORS 

The penetration and exit of flux computed from Eqs. 
(8-15) is visualized in Figs. 1-3 for isotropic strips and 
disks without volume pinning, using a flux-flow resistiv- 
ity pff — poB(r) with p = 140 (strip) or p = 70 (disk) 
in units where H c \ = a = pa = \dH a /dt\ = 1. Figure 1 
shows the field lines of B(x, y) in two pin- free strips with 
aspect ratios b/a — 2 and b/a = 0.3; Fig. 2 shows the 
surface screening currents in the same strips before flux 
has penetrated; and Fig. 3 plots some induction profiles 
in a strip and some hysteresis loops of the magnetization 
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FIG. 7. Magnetization curves of the same disk as in Fig. 
5 but for the Kim model, J C (B) = J c0 /(1 + 3|S|/B c i) for 
various pinning strengths J c o = 0, 0.25, 0.5, 1, 1.5, 2, 3, 4 in 
units H c i/a. Presentation as in Fig. 5. 

and of the induction in the center of a strip and disk. 

The profiles of the induction B y (r, y) taken along the 
midplane y — of the thick disk in Fig. 3 have a pro- 
nounced minimum near the edge r = a, which is the 
region where strong screening currents flow. Away from 
the edges, the current density J = V x B/|i is nearly 
zero; note the parallel field lines in Fig. 1. The quantity 
Jh = V x H(B) which enters the Lorentz force density 
Jh x B, is even exactly zero since we assume absence of 
pinning and the viscous drag force is small. Our finite 
flux-flow parameter po and finite ramp rate dH a /dt = ±1 
mean a dragging force which, similar to pinning, causes 
a weak hysteresis and a small remanent flux at H a = 0; 
this artefact is reduced by choosing a larger resistivity or 
a slower ramp rate. 

In Fig. 3 the induction Bq = B y (0,0) in the spec- 
imen center performs a hysteresis loop very similar to 
the magnetization loops M(H a ) shown in Figs. 3 and 
4. Both loops are symmetric, M(—H a ) = —M(H a ) and 
B (-H a ) = -B (H a ). The maximum of M(H a ) defines 
a field of first flux entry H en , which closely coincides with 
the field H' en at which B y (0, 0) starts to appear. The 
computed entry fields are well fitted by 

H s c f p /H cl = tanh y/0.36b/a , 

H^ sk /H cl = tanh ^Q.Q7b/a . (16) 

These formulae are good approximations for all aspect ra- 
tios < b/a < oo, see also the estimates of H Qn ps \fbfa 
for thin strips in Refs. jf]Jl(J]. 

The virgin curve of the irreversible M(H a ) of strips 
and disks at small H a coincides with the ideal Meissner 
straight line M = —H a /(\ — N) of the corresponding 
ellipsoid, Eqs. (4,7). When the increasing H a approaches 
H en , flux starts to penetrate into the corners in form of 
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FIG. 8. Magnetization curves as in Fig. 6 but for the Kim 
model Jc(-B) = J c o/(l + 3\B\/aJ c o) with J c o=const. for varios 
H c i — 0, 0.2, 0.35, 0.5, 0.7, 1 in units aJ c o- Also depicted are 
the pin-free magnetization (line with dots; M and H a here are 
in units H c i since J c o = 0) and the irreversible magnetization 
of the corresponding ellipsoid. 

almost straight flux lines (Fig. 1) and thus \M(H a )\ falls 
below the Meissner line. At H a — H cn flux penetrates 
and jumps to the center, and \M(H a ) \ starts to decrease. 
In decreasing H a , this barrier is absent. As soon as flux 
exit starts, all our calculated M(H a ) exhibit strong "nu- 
merical noise" , which reflects the instability of this state. 
Similar but weaker noise occurs at the onset of flux pen- 
etration. 

As can be seen in Fig. 4, above some field H rev , the 
magnetization curve M(H a ) becomes reversible and ex- 
actly coincides with the curve of the ellipsoid defined by 
Eqs. (3, 6, 7) (in the quasistatic limit with p^ 1 dH a /dt — > 
0). The irreversibility field H ICV is difficult to compute 
since it slightly depends on the choices of the flux-flow 
parameter po (or ramp rate) and of the numerical grid, 
and also on the model for Al(H a ;0). In the interval 
0.08 < b/a < 5 we find with relative error of 3%, 

H^ p /H cl = 0.65 + 0.12 In (b/a) , 

H^ k /H cl = 0.75 + 0.15 In (b/a) . (17) 

This fit obviously does not apply to very small 6/a«l 
(since H rev should exceed H cn > 0) nor to very large 
b/a ^S> 1 (where H ICV should be close to H c \). The lim- 
iting value of H lcv for thin films with b <C a is thus not 
yet known. 

Remarkably, the irreversible magnetization curves 
M(H a ) of pin-free strips and disks fall on top of each 
other if the strip is chosen twice as thick as the disk, 
(6/a) s trip ~ 2(6/a)disk- This striking coincidence holds 
for all aspect ratios < b/a < oo and can be seen from 
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FIG. 9. Same magnetization curves as in Fig. 5 but for a 
thinner disk with aspect ratio b/a = 0.125 for various degrees 
of pinning and constant H c \. 

each of Eqs. (7,16,17): The effective N [or virgin slope 
1/(1 — N)], the entry field H cn , and the reversibility field 
H ICV are nearly equal for strips and disks with half thick- 
ness, or for slabs and cylinders with half length. 

Another interesting feature of the pin-free magnetiza- 
tion loops is that the maximum of \M(H a )\ exceeds the 
maximum of the reversible curve (equal to H c i) when 
b/a < 0.8 for strips and b/a < 0.4 for disks, but at larger 
b/a it falls below H c \. The maximum magnetization may 
be estimated from the slope of the virgin curve 1 / (1 — N), 
Eq. (7), and from the field of first flux entry, Eq. (16). 

The formulae (7,16,17) are derived essentially from 
first principles, with no assumptions but the geometry 
and finite H c ±. They should be used to interpret exper- 
iments on superconductors with no or very weak vortex 
pinning. At present it is not clear how the presence of 
a microscopic Bean-Livingston barrier may modify these 
continuum theoretical results. 



VI. SUPERCONDUCTORS WITH PINNING 

Figures 5-8 show how the irreversible magnetization 
loops of disks with b/a — 0.25 (and in Fig. 9 for a thin- 
ner disk with b/a = 0.125) are modified when volume 
pinning is switched on. In Figs. 5, 6, and 9, pinning is de- 
scribed by the Bean model with constant critical current 
density J c , while in Figs. 7 and 8 the Kim model is used 
with an induction-dependent J C (B) — J c o/(l + 'S\B\/ Bk) 
with B K = (Fig. 8) or B K = n a aJ c0 /3 (Fig. 9). 

In Figs. 5, 7, and 9, H c i is held constant; with increas- 
ing J c or J c o (in natural units H c i/a) the magnetization 
loops are inflated nearly symmetrically about the pin- free 
loop or about the reversible curve (above H rev ), and the 
maximum of \M(H a )\ shifts to higher fields. Above H rev 
the width of the loop is nearly proportional to J c , as ex- 
pected from theories [[l5 16 which assumed H c i = 0, but 



at small fields the influence of finite H c i is clearly seen 
up to rather strong pinning. 

In Figs. 6 and 8, J c or J c q is held constant and H c i 
increased from zero (in natural units aJ c ). As expected, 
the influence of finite H c \ is most pronounced at small ap- 
plied fields H a , where it causes a peak in —M even in the 
Bean magnetization curves, which without consideration 
of H c i consist of two monotonic branches and a mono- 
tonic virgin curve. Within the Kim model, or with any 
decreasing J C {B) dependence, the magnetization loops 
exhibit a maximum even when H c \ = is assumed |l9[ ]. 
With increasing H c \ this maximum becomes sharper and 
shifts to larger fields, cf. Fig. 8. Comparing Figs. 5 and 9 
one sees that for superconductor disks with pinning and 
with H c i > 0, the peak in —M{H a ) becomes more pro- 
nounced and shifts towards smaller applied fields when 
the disk thickness is decreased. 
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